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Abstract 



We calculate the x-ray absorption near edge structure (XANES) modulating 
<^> the Ba 4d-photoabsorption cross section - the giant dipole resonance - centered 

I around 110 eV photon energy for several models of spherical and non-spherical 

Ba@C n metallofullerene cage systems. In the cases considered, the XANES 
interference patterns provide clear structural "fingerprints", distinguishing 
between center versus off-center Ba position and spherical versus deformed 
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C n shell. 
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Cage molecules are well known for spectacular interference effects, causing very pro- 
nounced structure in photoabsorption cross sections. Nearly ideal cages should be provided 
by fullerene molecules C n , n — 60, 70, . . ., in gas phase or in solid form. Among various 
forms of doped fullerenens, there is the possibility to put metal atoms inside C„, endohedral 
metallofullerenes - M@C„ (!]]. Recently there has been a wealth of reports on the preparation 
and characterization of endohedral metallofullerenes M@C n as well as calculations of 

structure, affinity and ionization energies PJ|,0]. 

The presence of a heavy metal atom inside a cage of ligands provides a very interest- 
ing many-electron system: the collective dynamics of the coupled many-electrons systems 
may give rise to giant dipole resonances, and single-particle multiple scattering will lead to 
standing wave patterns and resonances inside the cage [@||. As a result, the broad giant 
dipole resonance will be modulated by x-ray absorption near edge structure (XANES) as 
illustrated by the Ba 4d-photoionization cross section in solid BaF 2 and YBaCuO (see e.g. 
ref. gg). 

In this Letter we present results from an investigation of the collective and single-electron 
response in photoemission from the 4d core level of a Ba atom placed inside a model fullerene 
cage in cases of spherical and non-spherical symmetry. We focus attention on two situations: 
(i) a Ba atom displaced from the center of a spherical Ceo model cage, Ba@ Ceo (Coot,), and (ii) 
a Ba atom at the center of a non-spherical cage consisting of a short tube with hemispherical 
endcaps, to be referred to as Ba@C 90 (D oo /j). We find that the multiple scattering structure 
of the photoabsorption cross section is a sensitive probe of the position of the emitter atom 
and the shape of the cage. We propose that the oscillations of the 4d-f giant dipole resonance 
can provide a chacteristic "fingerprint" of the system - distinguishing almost by inspection 
between center vs off-center position of the Ba emitter atom as well as spherical vs non- 
spherical shape of the cage. 

The Ba@C n structure is defined by a point charge Ze a = 56 for the Ba nucleus and 
a surface effective nuclear charge Z she \\ for the C n shell ||. In this way the nuclear point 
charges of the real C n structure is averaged to a homogeneous surface charge distribution 
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corresponding to the number of carbon valence electrons (4 per C-atom). This is in fact 
an effective 2-dimensional jellium model. We then fill in the appropriate number iV = 
Ze a + Z S h e n of electrons for a neutral Ba@C n cluster and calculate the electronic structure 
selfconsistently within the local-density approximation (LDA) (the structure and dynamics 
of spherical La@C 60 was studied in the same way by Wastberg and Wendin ||). 

To find the electronic structure of a non-spherical cluster we expand the initial bound 
state wave functions ^(r) = (1/r) J2l UL^Y^ft) and the final state Green functions 
G(r,r',E) = J2L,L'YL(^)GL,L'(r,r')Yl,(Q'), L = (l,m), in spherical harmonics around the 
center of the heavy central atom JTT]JT^] . In such a one-center expansion, the molecular 
problem reduces to an atomic problem, but with coupled angular-momentum channels. The 
advantage is that the Coulomb interaction can be expressed in terms of spherical harmonics, 
which allows straigthforward calculation of the dielectric response function. Hence electron- 
electron correlation and many-particle dynamics in the cluster can be studied using well 
established many-body techniques for atomic systems. 

The resulting coupled radial Schrodinger equations were solved by procedures described 
in Ref. flT2|. From the spherical expansion of initial state wave functions and final state 



Green functions |T3] we then construct the dielectric susceptibility xo( r , r ';^) and calculate 
the induced charge 8n(r; u) in a self-consistent manner within the TDLDA (time depen- 



dent local density approximation) |14|. We finally obtain the photoabsorption cross section 



o~(uj) ~ ulm J z5n(r;uj)dr where z = rcos6 is the the electric dipole operator [15j with the 



electric field along the z-axis. In this work the polarization vector is always chosen along 
the rotational symmetry axis |16] . 



Figure [I] shows the electronically selfconsistent (SCF) (fixed geometry) ground state 
charge densities for the three different model systems we consider. Figure |l|(a) shows the 
electronic density of spherically symmetric Ba@C 60 , Fig.|l](b) shows the case of Ba@C 60 (C ool ,) 
with Ba displaced from the center of a spherically symmetric Ceo cage, and Fig.|l|(c) shows 
the symmetric case of Ba@Cgo(D 00 / l ), where the cage is a short cylinder with hemispherical 
Ceo endcaps. The radius of the Ba 5s/5p charge density is about 2.5 a.u., corresponding 



to the radius of the well isolated central charge density of the Ba 2+ ion core (black spot) 
in Figs.(la)-(lc). In contrast, the Ba 6s-radius is about 6 a.u., causing large 6s-amplitude 
beween Ba the C n shell. The calculated low charge density in this bond region suggests 
that the two Ba 6s valence electrons have been transferred to the C n shell, leaving a doubly 
ionized central Ba 2+ ion, (Xe 4c? 10 5s 2 5p 6 configuration), with some weak effects of bond 
formation in the displaced case in Fig.^(b). In the case of Ba@Cgo in Fig.^(c), the strongly 
ionic character is particularly evident. 

The geometries in Fig{l] do not necessarily represent equlibrium positions. In fact, we 
have calculated the variation of the total energy for Ba@Ceo as a function of displacement of 
the Ba atom. Technically, the shell is displaced, and we study the variation of -E to t(Ba@C6o) — 
-Etot(C6o) to minimize systematic errors from the finite /-expansion. As a result we find that 
there is a small lowering of E tot when the Ba atom is displaced from the center of Ceo- Since 
we have not allowed any freedom in the shape and structure of the Cqq shell, we cannot 
draw any firm conclusions. However, our results suggest that here is a tendency for the Ba 
atom to have an off-center equilibrium position, in agreement with calculations for La@Cs2 
|J, Ce@C 82 |7j and Ca@C 60 ||. It might even be appropriate to regard the metal atom as 
"adsorbed" on the inside of the fullerene cage. 

Figure |2| shows a central result of this Letter, namely the photoabsorption (total pho- 
toionization) cross sections <y{uo) of the Ba@C n clusters described in Fig. 1. The electric 
field vector is along the rotational symmetry axis in displaced Ba@C6o(Coo,;) and along the 
cylinder (long) axis in elongated Ba@C 90 (D oo / l ). 



From our previous work [§,|TlJ] we know that the environment does not modify the overall 
strength and shape of the Ba 4d — ef (4d-"4f ) giant dipole resonance: the effect is to induce 
a structure which oscillates around the atomic "background". The three curves in Fig.|2| 
nevertheless present three easily distinguishable patterns - or "fingerprints" - which may 
serve as a guide to the structure of the local environment: 

(i) The spherical Ba@C 60 case (Fig.|l](a)) is simple: the pronounced XANES oscilla- 
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tions present a clear signature of the unique distance of a single coordination shell. These 
oscillations must be associated with additional nodes of the continuum e/-wave function 
moving through the Ceo shell into a radial standing wave in the inner region outside the Ba 
core. This f-wave function must be of a molecular, multiple scattering, kind connected with 
large-amplitude shape resonances, and cannot be described in terms of single scattering. 

(ii) The non-spherical, displaced, Ba@C6o(Coo„) case (Fig.|l|(b)) lacks pronounced oscil- 
lations and, at first sight, even characteristic structure seems to be absent. However, the 
shape of the 4d cross section is quite different from that of a free Ba atom superimposed on 
a Ceo background [Q (which corresponds to averaging out the oscillations in Ba@Ceo in case 
(i) above). With a displaced Ba atom there are now obviously two characteristic distances 
along the electric field vector: a short "bond" distance of about 5 a.u. and a long distance 
of about 8.4 a.u. A closer inspection of Fig.|^ suggests that the XANES peaks in case (i) 
have been shifted towards higher energies, perhaps due to the a dominating influence of the 
shortest Ba — C 60 "bond" distance. 

(iii) The non-spherical, but inversion symmetric, Ba@Cgo(D 00 / l ) case (Fig.|I|(c)) lacks 
pronounced XANES oscillations but nevertheless shows considerable characteristic structure. 
There is no single unique coordination distance, but - as will be demonstrated below - one 
can observe structure due to standing waves along the cylinder (long) axis. 

To illuminate the electric-field polarization dependence of the XANES we compare the 
cylindrical, cigar-shaped Ba@Cg (D oo / l ) with the corresponding inscribed (short axis) and 
circumscribed (long axis) spherical systems. The inscribed system we have already discussed 
- Ba@Cgo- The circumscribed model system is provided by a spherically symmetric Ba@Ci3o 
cluster (522 electrons in the isolated jellium shell). Ba@Ci3o is designed to have a unique 
coordination length equal to the long axis of Ba@Cgo - the question is whether this similarity 
will show up in the XANES. 

In Fig. |3], the similarity between Ba@Ci3o and Ba@Cgo is striking: the peak portion 
of the Ba@Cgo cross section looks very similar to that of Ba@Ci3o with the oscillations 
averaged out. This suggests that the flat-topped giant 4d-dipole resonance of Ba@Cgo does 
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reflect the narrow distribution of coordination distances defined by the long axis. Moreover, 
since Ba@Ceo corresponds to the short axis and Ba@Ci3o to the long axis of Ba@Cgo, Fig.|3| 
suggests that there will be constructive interference at the Ba nucleus (4d-shell) around 110 
and 140 eV photon energy and destructive interference around 120 -130 eV, as indeed seems 
to be the case in the Ba@C 90 cross section. 

Figure ||| also allows comparison of the XANES of the spherical Ba@Ceo and Ba@Ci3o 
systems. The immediate impression of the oscillating structures is that they arise from 
standing waves inside the shells, and that the wavelengths must corespond to the radii of 
the shells. Plotted on a final state momentum scale k = y/e = \/u> + E^, in both cases 
the peaks appear as roughly equally spaced above u — 110 eV. In the simplest of models 
we may try to associate the peaks with high-lying resonances in a radial quantum well, 
kR = nir + const, R = ir/5k. In this way R(Ba@C 60 ) « 7-9 a.u. and R(Ba@Ci 30 ) « 9-11 
a.u., in reasonable agreement with the given radii of the shells (6.7 and 10 a.u, resp), but 
with considerable variations from peak to peak. Therefore, even if the XANES oscillations 
look nice and regular, they cannot be used for precise structure determination in the spirit 
of high-energy EXAFS (extended X-ray absorption fine structure). 

We are not aware of any experimental results for endohedral metallofullerene systems 
that can be compared with the present theoretical photoabsorption cross sections. However, 
it is a fact that the 4d-giant dipole resonance in the high-T c superconductor YBaCuO (Ba 
in a cage of 6- fold coordinated O; Ba-0 distance about 5.1 a.u.) looks very similar to the 
case of spherically symmetric Ba@Ceo || shown in Fig.|l|. To be able to observe such strong 
oscillations in metallofullerens, according to the present investigation the metal atom must 
sit at the center of a spherical cage (unique coordination distance). Since this is not likely 
to occur for Ba@Cgo, because the cage is too large and the Ba atom will sit at an off-center 
position, it might be necessary to laser evaporate the C n fullerene cage and shrink it down 
to, say, Ba@C44 (in analogy with La@C44 Q). In this case the cage radius (R ~ 5.7 a.u.) 
will correspond more closely to a natural chemical bonding distance - possibly this is the 
case that will compare most directly with the XANES-modulated 4d-giant dipole resonance 
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in YBaCuO. 

Finally, in Fig.|] we present partial photoionization cross sections of displaced 
Ba@C6o(C 001) ). As expected, the 4d-f giant dipole resonance shows upp in all of the Ad, 5s, 5p 
and valence level photoemission cross sections with their usual characteristic resonance pro- 
files: n = 4 subshells roughly follow the 4d-f resonance, while n> 5 subhells (which are 
basically outside the 4d-f dipole) show their resonance enhancement on the rising edge of 
the 4d-f giant dipole resonance. The weak hump at around 125 eV, discussed above, can 
be seen in all of the emission channels. The valence emission shows resonance enhancement 
around 100-110 eV at the same position as the emission from the 5s and 5p orbitals, and 
then goes over into "background" emission from the Cgo shell. This indicates significant hy- 
bridization of the shell with Ba 5d orbitals, in line with the recent photoemisson results 
for solid Ba 6 C 60 by Knupfer et al. |T7|JT6|| and with calculations for Ba 6 C 60 UP9pO| . 



In fact, our complete results show that the two highest ocupied molecular orbitals 
(HOMO) have 4f-character on Ba while the next lower orbitals have 5d- and 6s/6p-character. 
This follows in two ways: (i) direct inspection (or partial-wave analysis) shows that the va- 
lence orbital character on Ba is mainly that of (distorted) Ba 4f, 5d and 6s/6p; (ii) the 
partial cross sections show resonance enhancement of the kind demonstrated in Fig.^ - the 
one for the 4f-hybridized orbitals looks similar to the 4d-cross section, while the ones for the 
5d- and 6s/6p-hybridized orbitals look similar to the 5s/5p cross sections. In this way we 
have a very useful connection between the character of the initial orbital of the photoemitted 
electron and the frequency dependence of the photoelectron intensity. 

In conclusion, we have applied a TDLDA one-center expansion to a model of a non- 
spherical metallofullerene cage system - Ba@C n - in order to study the combined effects of 
dynamic screening, multiple scattering and non-spherical local environment on the collective 
Ba 4d-giant dipole resonance centered around 110 eV photon energy. For the geometries 
selected in this work, the XANES show distinctly different patterns, distinguishing almost 
by inspection between center versus off-center position of the Ba emitter atom as well as 
spherical versus non-spherical shape of the cage. Our results directly concern photoabsorp- 



tion and photemission from oriented samples - e.g. solids or adsorbates - with polarization 
vectors along the high-symmetry rotation axis. However, our analysis suggests that sig- 
nificant parts of the corresponding XANES will be recognizable also in gas phase spectra. 
Whether these XANES signatures are robust enough to serve as unambiguous " fingerprints" 
in general situations can only be decided by comparison with future experimental results. 
It would be useful to have photoionization cross sections for Xe@C n f21J and surrounding 
elements, e.g. Sn@C„-Ce@C n , representing different cases of bonding and charge transfer, 
shapes of 4d-giant dipole resonances, and many-electron effects. 
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FIGURES 

FIG. 1. Self-consistent charge density plots for Ba@C n metallofullerene model clusters (2D 
jellium model for the C n nuclear charge; number of electrons chosen to give closed-shell electronic 
structure): (a) Ba@C6o: spherically symmetric Ceo cage (R=6.7 a.u.; 56+230 electrons) with Ba 
at the center; (b) Ba@C6o(Cocn;): spherically symmetric Ceo cage (R=6.7 a.u.; 56+230 electrons) 
with Ba displaced 1.725 a.u. from the center); (c) Ba@Cgo(D 00 / l ): Ba at the inversion center of 
C n cylinder with hemispherical Geo end caps (-R cy i = 6.7 a.u., total length=10.0 a.u.; 56+360 
electrons). 

FIG. 2. Photoabsorption cross sections of the Ba@C n clusters defined in Fig. 1: Spherically 
symmetric Ba@Cgo (E^ =90AeV); displaced Ba, Ba@C6o(Coo»;) {Em =89.6eF); elongated C n , 
Ba@Cgo(D 00 / l ) (E^d =91.6eV). For comparison, = 93.7eU in a free LDA Ba atom, and 



E Ad =93.3eV in selfconsistent (SCF) O-Ba-0 (the result in Ref. Q is not SCF). All calculated 
energies in this work are non-relativistic. 

FIG. 3. Photoabsorption cross sections of Ba@C n clusters with an inversion center: Spherical 
Ba@C6Cb non-spherical Ba@Cgo(D 00 /j), and spherical Ba@Ci3o {E±d = 94.3eF). 

FIG. 4. Partial photoionization cross sections of displaced Ba@C6o(Cooi,) (see also Fig. 1(b) and 
Fig. 2). The 4d cross section is the sum of 4da, 4dir, and Ad5 components with energy splitting 
0.016 eV (« 0.055 eV in elongated Ba@Cgo). The 5p cross section (E§ p = 16.4eV) is the sum 
of bpa and bpir components with energy splitting ~ 0.116 eV. The 5sa level has binding energy 
E$ s = 28.2eV. The valence cross section represents the integrated emission from the occupied 
levels in the 3.14(HOMO)-8.60 eV range below the vacuum level. 
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Partial Cross Sections Ba@C 
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